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. Manifestations of pronounced shell effects are discovered when adding 

^ . nonaxial octupole deformations to a harmonic oscillator model. The degen- 

. eracies of the quantum spectra are in a good agreement with the corresponding 
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main periodic orbits and winding number ratios which are found by classical 
analysis. 
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Q>^ . many investigations into the contribution of higher multipoles on the formation of shell 



The remarkable regularity of rotational spectra of superdeformed nuclei has prompted 



structure [Q. Related questions have arisen for other mesoscopic systems. In particular, it 
appears that octupole deformation has the same importance for two-dimensional systems 
like quantum dots and surface clusters |^ as for three-dimensional systems like nuclei ^ and 
^ , metallic clusters [Q. Semi-classical analysis based on periodic orbit theory |^-|^ provides 
substantial insight into the role of the octupole deformation for axially symmetric systems 
||8HlO|. Axially symmetric octupole deformation which leads to soft chaos in the classical 
case produces short periodic orbits at particular parameter strengths, and, correspondingly, 
pronounced shell effects arise in the quantum spectrum @J^. 

Conservation of angular momentum may increase the regular region for a non-integrable 
problem with axial symmetry (see for example The situation becomes more compli- 



cated for nonaxial systems with three degrees of freedom, since angular momentum is no 
longer a constant of motion and the classical dynamics may lead to a stronger degree of 
chaos |T^. Inclusion of exotic, i.e. nonaxial, octupole deformations renders the finding of 
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pronounced shell effects rather difficult. Results based on the term Isii, which was sug- 
gested by Mottelson [^] and studied in [Q, have been questioned in |jl5[. Other attempts 



which incorporate nonaxial octupole deformation start with axially symmetric potentials; 
they have found indications of shell effects using deformations mainly with /i = 0, 2 
|]T^,|TB|. The increasing accuracy of measurements of nuclear spectra, due to the new gener- 
ation of detectors, gives substantial indications of strong octupole correlations ||3|,p!7|. This 



calls for a thorough analysis of nonaxial octupole deformations. A similar question about 
shell effects in electronic structures and their connection to octupole deformations arises 



from the ab initio calculations of the melting transitions in small alkali clusters [0. In 
this Letter we demonstrate the existence of strong shell effects which arise in the tri-axial 
harmonic oscillator combined with nonaxial octupole deformations. The model may serve 
as a simple and transparent study of the effective mean field for mesoscopic systems like 
nuclei and metallic clusters. In addition, features characteristic of realistic potentials, i.e. a 
coexistence of regular and chaotic dynamics and the consequences for quantum mechanics, 
are addressed. 

Similar to the procedure pursued in ^ we are guided by the study of the classical motion 
in obtaining the quantum mechanical results. The single-particle Hamiltonian considered 

reads 



H = ^ + 
2m 



fj. . 
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where the Ya^j are the usual spherical harmonics. To maintain time reversal invariance, only 
the combinations Y^^ ± Y3_^ are considered, with a factor i where appropriate. We take into 
account only one of the deformation /i = 0, 1, 2, 3 at a time. For convenience, we express all 
quantities below in units of uj = uj^, i-e. a = 1. If any of the parameters As^ are nonzero, we 
are faced with a non-integrable system. In fact, the problem gives rise to chaotic motion even 
at relatively small values of the octupole parameters. The parameters have to be limited 
by their respective critical values Ag^** at which the potential no longer binds. The critical 
values depend on the parameters a, b, c which can be expressed through standard quadrupole 



deformations and, if more than one Aa^ is considered, critical surfaces are obtained. 
It is obvious that a search for shell structure for the corresponding quantum mechanical 
problem becomes meaningless above the critical values as then the quantum mechanical 
spectrum obtained by matrix diagonalisation does not relate to the corresponding classical 
Hamiltonian. 

The quadrupole shapes as determined by the parameters a, b, c are illustrated by the 
hexagonal figure given in Fig.l. The three axes denoting axially symmetric prolate (oblate) 
shapes differ by an appropriate permutation of the coordinates x — y — z. For physical con- 
sideration it is therefore sufficient to consider just one sector if only quadrupole deformation 
is being studied . However, the addition of an octupole term defines an orientation, since 
it refers to a specific z-axis. Therefore, when adding, say, the term t'^X^qY^q, the physical 
situation is different for two oblate cases, for example, obl/x and obl/z in Fig.l. Note that 
the latter case would preserve axial symmetry, thus making it effectively a two degrees of 
freedom system. But the former case is now a genuine three degrees of freedom system, 
since the symmetry axis of the quadrupole shape (the x-axis) is different from the symmetry 
axis of r^Asolso (the 2;-axis). 

In the same vein, the addition of a term r^As^Ysi^ to any of the axially symmetric 
quadrupole shapes, gives rise to a three degrees of freedom system for /i 7^ 0. 

The effect of r'^Y^Q upon the axes pro/ z and obl/z has been dealt with in We here 
begin with the combination r'^iY^-^ — Y^^^). It appears from numerical analysis of the classical 
equations of motion (see below) that the most regular motion occurs in the vicinity of the 
pro/x axis. The procedure to approximate the non-integrable classical system is the 'removal 
of resonances' method [1^. To lowest order it consists of averaging the Hamilton function 
over the fastest angle of the unperturbed motion (all As^ = 0) after rewriting the momenta 
and coordinates in terms of action-angle variables, viz. 



= Pi = V^JimUJi COS 6i, 

6i = uJit, i = x,y,z. 
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On the axis pro/x one would therefore average over 9z or 9y and above it over since there 
ujz > ujy > ujx {a > b > c). Based on the observation that r^iY^^ — Y^^^) is proportional to 
a;(x^ — 3y'^)/r, we expect the z-motion to be weakly affected by this term. In the vicinity of 
the pro/x axis an averaging over 6z is thus indicated. Moreover, for the reasons just given, 
an averaging over 6z, which yields an effective potential in the x — y-coordinates, is expected 
to make little difference from an effective potential obtained by simply setting z = 0. This 
expectation is convincingly confirmed by numerical tests as long as Uz > ujy. In other words, 
for ujz > ujy > ujx, the motion effectively decouples into an unperturbed motion in the z- 
coordinate (governed by the potential ■mz^uj'^/{2c^)) and the two degrees of freedom motion 
in the x — y-plane. Averaging now over the fast angle 6y yields the unperturbed motion in 
the ^/-coordinate (governed by m?/^u;^/(26^)) and the effective potential for x which reads 



U,s{x) = ^[x^ + 
A33^(2x2K(-||) - 3<,22Fi(i, f , 2; -||))] (2) 

where K and 2-^1 denote the first elliptic integral and the hypergeometric function, respec- 
tively, and = 2Ey/ {mojy) where Uy = u/b. The approximation used here assumes that 
Ey, the energy residing in the y-motion, is constant (and therefore also E^); note that the 
effective potential Ucs{x) depends on Ey. A numerical comparison between the true three 
dimensional motion and the approximate decoupled motion nicely confirms the validity of 
the approximation for A33 < Acrit/2. The crucial test which is relevant for the correspond- 
ing quantum mechanical case is the comparison of the winding number ratios ul^/uy and 
ul^/uz- Moreover, these ratios are virtually independent of Ey for Ey less than 60% of its 
maximal value -Etot- As a result, we may evaluate analytically by choosing i^^y = in 
Eq.(^ and obtain 



yi - 1 1 ... 

V1-A33 + V/1 + A33^&^C 

By appropriately tuning the parameters b, c and A33, we can obtain simple ratios for the 
winding numbers. They determine short periodic orbits which are expected to occupy a 



major part of a phase space. 

For a two-dimensional problem the Poincare surface of section can be used for the esti- 
mation of the percentage of phase space occupied by periodic orbits. In the present situation, 
the five dimensional phase space renders an understanding of the underlying structure rather 
difficult. The technique utilized here is essentially a frequency analysis. It is in principle 
impossible to determine whether a trajectory is quasi-periodic or chaotic merely by looking 
at the frequency spectra [Q. A pragmatic approach is adopted here in that an initial con- 



dition is deemed to yield a chaotic orbit if the associated frequency spectra have sufficiently 
many discernible peaks. We call an orbit chaotic if any one of the frequency spectra has 
more than six peaks with intensities greater than 1% of the maximum intensity. These 
arbitrary choices proved satisfactory for our purpose. If the orbit is quasi-periodic, then 
the most significant frequency peaks are compared, and the approximate winding number 
and period obtained. Repeating this procedure many times using different initial condi- 
tions in phase space yields a Monte-Carlo type estimate of the portions of phase space 
characterized by the various different frequency ratios. If a particular simple ratio domi- 
nates, as exemplified below, then we expect specific signatures in the quantum spectrum 
as shell structure. Examples are illustrated in Figs. 2-4. We have chosen the numbers 
a : 6 : c = 1 : 0.5577 : 0.5577, 1 : 0.3718 : 0.5577 to obtain the ratios 1 : 2 : 2, 1 : 3 : 2 from 
Eq.@, respectively. This ratio is sufficiently simple to make an easy comparison with the 
quantum results. 

The adiabatic approach predicts orbits with the frequency ratios 1:2:2 and 1 : 3 : 2 at 
A33 = 0.5Xcrit- The corresponding spectra are displayed in Fig.??. The classical frequency 
analysis (CFA) of the exact orbits shows a peak at A33 ~ O.55A33** for the ratios 1:2: 2,1 : 
3 : 2, and the quantum shell structure occurring at A33 ~ O.5A33** has the correct degeneracy 
pattern for about the first hundred levels. 

As a quantitative measure for shell structure we use the Strutinsky-type analysis intro- 
duced in §. From the quantity AE{\,N) = 5E{\,N + 1) + SE{\,N - 1) - 25E{\,N), 
where 6E is the fiuctuating part of the total energy, we obtain the precise location of the 



magic numbers (see Fig.4). Similarly, the whole discussion can also be applied to the axis 
pro/y in Fig.l by using the combination r^(Y33 + Ys^s) ~ y{y^ — 3a;^)/r instead. It appears 
that, due to the weak z-dependence of the combinations r^(Y33 ± Ys-s), mainly the x — y 
profile of the unperturbed harmonic oscillator is important. The stronger the deformation 
in the x — ?/-plane, i.e. the further away from the 2;-axial symmetry line {pro/ z and obl/z), 
the better the adiabatic approximation and shell structure becomes. In contrast, along this 
horizontal line either l3±3 combination acts upon a circular potential in the x — y-plane and 
quickly introduces chaos. 

The CFA, applied to the combination r'^{Y^i — Ya-i) ~ x{4:z'^ — — y^), reveals that 
the y motion is weakly affected for A31 < O.GAg'J** in the vicinity of the pro/x axis. Applying 
the analysis described above, we found that this term leads to shell effects similar to those 
of r^(y33 — 13-3). In fact, Eq.(^ holds for this combination as well. The plus combination 
^"^(^31 + ^3-1) is simply the minus combination under the interchange of x and y, and thus 
will produce the same effects near the region pro/y. It is important to note that the addition 
of a term X-^^r'^Y^Q to either situation, pro/x or pro/y, leads to an onset of chaos for rather 
small values of A30, and accordingly the quantum spectrum does not exhibit shell structure. 

The two combinations r'^iY^^ -^Y^^^) produce the same potential shape for the spherical 
case [|l6l. However, according to the analysis above, the plus and minus combinations 
have different effects for different sectors of the hexagonal figure (Fig.l). The CFA shows 
that the adding of the octupole term A32r^(F32 + I3-2) gives rise to chaotic motion for 
comparatively small coupling values. In contrast, the term A32r^(y32 — 13-2) has less impact 
on the unperturbed motion, chaotic motion only becoming discernible for A32 > O.SAcrit- 
Since (F32 — 1^3-2) is symmetric with respect to an interchange of x and y, this result 
applies to the region above and below the axis pro/z including the axes ohl/x and ohl/y. 
The quantum mechanical results are in accordance with the classical findings: the plain 
quadrupole spectrum changes weakly over a considerable range of A32 when adding the term 
'^^{Yz2 — ^3-2) while the order soon decays when the plus combination is switched on. 

The cases considered here represent novel examples of three-dimensional non-integrable 



systems, which can be well-approximated by integrable ones. The results of the current 
literature are limited primarily to axially symmetric non-integrable systems ||7|,p!2|. Special 
values of parameters, found with the 'removal of resonances' method, produce potentials 
conducive to regular classical motion in much of phase space. The various octupole combi- 
nations may have different effects on the generation of shell structure, depending on where 
the unperturbed potential lies in the hexagonal figure (Fig.l). The effect of r'^Y^o upon the 
axis pro/ z is similar to that of r'^{Y^-i ± I3-3) and r'^iY^i ± F3-1) upon the axis pro/x or 
pro/y. In contrast, the terms r'^{Y^2 ± ^3-2) do not support shell structure. In this context 
we mention that the special combination r'^XiY^Q + 8(132 + 13-2)), when added to pro/ z, 
is, after suitable permutation of the coordinates, identically equivalent to the adding of 
r^A33(y33 — 1^3-3) to pro/x. Thus, we have generalized our previous result [Q to the domain 
of triaxiality in that the combination of quadrupole and nonaxial octupole deformations has 
been shown to lead to shell effects equivalent to those from a plain quadrupole deformed 
potential, at least for the first one hundred levels. More comprehensive details of the classi- 
cal frequency analysis, the adiabatic approach and the quantum mechanical analysis will be 
presented in a forthcoming paper. Finally, we note that negative parity states observed in 
rare-earth nuclei with neutron number ~ 92, which become yrast at high spins, need an 



unexpectedly strong degree of triaxiality (see [^]) when described in terms of quadrupole 
and hexadecapole deformations only. We suggest that inclusion of an octupole deformation 
^"^(^33 ~ ^3-3) or ^ banana-type octupole deformation r^(y3i — Y^_i) which effectively gives 
rise to shell effects of a tri-axial oscillator, could yield a more natural explanation of these 
phenomena. 
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Figure Captions 

Fig.l Shapes in the a, b, c plane. Spherical symmetry (a = 6 = c) is at the center while 
axially symmetric prolate and oblate shapes are obtained along the various axes. A genuine 
tri-axial quadrupole deformation {a ^ b ^ c) occurs between the axes. 

Fig.2 Two quantum spectra: (a) for the parameters b — c — 0.5577; (b) for the param- 
eters b — 0.3718, c = 0.5577. Energies are given in units of ku;. 

Fig.3 Classical phase space occupation for the frequency ratios 1.1:3:2 and 1:3:2. The 
1:3:2 peak at A33/A33** f» 0.55 leads to the shell structure displayed in Fig.2b. 

Fig.4 Magic numbers calculated at values of parameters as for Fig.2. The magic numbers 
corresponding to a pure quadrupole deformed harmonic oscillator are indicated in (a) for 
the ratio 1:2:2 and in (b) for the ratio 1:3:2. 
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